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Abstract 

Recently it is shown that there are three families of stochastic one-dimensional 
non-equilibrium lattice models for which the single-shock measures form 
an invariant subspace of the states of these models. Here, both the station- 
ary states and dynamics of single-shocks on a one-dimensional lattice are 
studied. This is done for both an infinite lattice and a finite lattice with 
boundaries. It is seen that these models possess both static and dynamical 
phase transitions. The static phase transition is the well known low-high 
density phase transition For the ASEP. The BCRW, and AKGP models 
also show the same phase transition. Double-shocks on a one-dimensional 
lattice are also investigated. It is shown that at the stationary state the 
contribution of double-shocks with higher width becomes small, and the 
main contribution comes from thin double-shocks. 
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1 Introduction 



Reaction-diffusion systems, is a well-studied area. People have studied reaction- 
diffusion systems, using analytical techniques, approximation methods, and 
simulations. A large fraction of exact results belong to low-dimensional (spe- 
cially one-dimensional) systems, where solving low-dimensional systems shoiild 
in principle be casi(T. D(^spite tlicnr simplicity, these systems exhibit a rich and 
rather non-trivial dynamical and stationary behavior. Studies on the models far 
from equilibrium have shown that there is a remarkably rich variety of critical 
phenomena [1]. 

Shocks in one-dimensional reaction-diffusion models have been absorbed 
many interest recently [2-12]. There are some exact results about shocks in 
one-dimensional reaction-diffusion models together with simulations, numeric 
results [7] and also mean field results [3]. Formation of localized shocks in 
one-dimensional driven diffusive systems with spacially homogeneous creation 
and annihilation of particles has been studied in [13]. Recently, in [5], the 
families of models with travelling wave solutions on a finite lattice have been 
presented. These models arc the Asymmetric Simple Exclusion Process (ASEP), 
the Branching- Coalescing Random Walk (BCRW) and the Asymmetric Kawasaki- 
Glauber process (AKGP). In all of these cases the time evolution of the shock 
measure is equivalent to that of a random walker on a lattice with L sites 
with homogeneous hopping rates in the bulk and special reflection rates at the 
boundary. 

Shocks have been studied at both the macroscopic and the microscopic levels 
and there are some efforts on addressing the question that how these macroscopic 
shocks originate from the microscopic dynamics [8]. Hydrodynamic limits are 
also investigated. 

Among the important aspects of reaction-diffusion systems, is the phase 
structure of the system. The static phase structure concerns with the time- 
independent profiles of the system, while the dynamical phase structure concerns 
with the evolution of the system, specially its relaxation behavior. In [16- 
19], the phase structure of some classes of single- or multiple-species reaction- 
diffusion systems have been investigated. These investigations were based on 
the one-point functions of the systems. 

Here we study both stationary and also dynamics of the single-shocks on a 
one-dimensional lattice. This is done for both an infinite lattice and a finite 
lattice with boundaries. In the stationary state, the system can be found in the 
low-density or in the high-density phase. This phase transition is a well known 
first order phase transition in familiar ASEP, and also its extensions [20]. The 
BCRW, and AKGP models show the same phase transitions. We also investigate 
the dynamical phase transitions of the models. It is seen that ASEP has no 
dynamical phase transition, but both of the models BCRW, and AKGP have 
three phases, and system may show dynamical phase transitions. Double shocks 
on a one-dimensional lattice have been also investigated, and its the stationary 
behavior have been studied. It is shown that, in the thermodynamic limit, 
contribution of double shocks with higher width become vanishingly small. 
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2 Fixing the notations 

Consider a one-dimensional lattice, each point of which is either empty or con- 
tains one particle. Let the lattice have L sites. An empty state is denoted by 
|0) and an occupied state is denoted by 

io,.(;), ,1, (1, 

If the probability that the site i is occupied is pi then the state of that is 

represented by ^"'^ Z'*^ • The observables of a reaction-diffusion system are the 

operators TV", where i with 1 < i < L denotes the site number, and a = 0, 1 
denotes the hole or the particle: Nf is the hole (vacancy) number operator at 
site J, and is the particle number operator at site i. One has obviously the 
constraint 

= 1, (2) 

where (s| is a covector the components of which (sq's) are all equal to one. The 
constraint 101, simply says that each site is either occupied by one particle or 
empty. A representation for these observables is 

iVf l0- - -0l (8)jV° ® 1 ^ ■ ■ ■ ^ 1 , (3) 

i— 1 L — i 

where N"" is a diagonal 2x2 matrix the only nonzero element of which is the 
a'th diagonal element, and the operators 1 in the above expression are also 2x2 
matrices. The state of the system is characterized by a vector 

|P) e y ■ ■ ■ V , (4) 

L 

where V is a 2-dimensional vector space. All the elements of the vector |P) are 
nonnegative, and 

(S|P> = 1. (5) 

Here (§| is the tensor-product of L covectors (s|. The evolution of the state of 
the system is given by 

|P)=7^|P), (6) 

where the Hamiltonian Ti is stochastic, by which it is meant that its nondiagonal 
elements are nonnegative and 

(§1 7i = 0. (7) 
Two conventions are used to write the master equation, 

a) |P)=-H'|P) ^, |P)(t) =exp(-m') |P)(0) 

b) |P> = n |P> ^, |P)(t) = exp [m) |P>(0) 
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These two conventions are related to each other simply through H' = —H. In 
the a) convention nondiagonal elements of Ti.' are negative of reaction rates, 
hence nonpositive, and its diagonal elements are nonnegative. In the b) conven- 
tion nondiagonal elements of Ti are reaction rates, hence nonnegative, and its 
diagonal elements are nonpositive. In the a) convention the real parts of the 
eigenvalues of Ti' are nonnegative, and the eigenvalue with minimum nonzero 
real part corresponds to the relaxation time. In the b) convention the real 
parts of the eigenvalues of H are nonpositive, and the eigenvalue with maximum 
nonzero real part corresponds to the relaxation time, and finally the state vector 
|P) is the same in two conventions, and in both conventions the elements of |P) 
are nonnegative. Through out this paper we use the b) convention. 

The interaction is nearest-neighbor, if the Hamiltonian is of the form 

L-l 

'H = ^ i^^j+i, (8) 

i=l 

where 

:^ 1 (8) ■ ■ ■ ^ 1 ®H® 1 (9) 

(It has been assumed that the sites of the system are identical, that is, the 
system is translation-invariant. Otherwise H in the right-hand side of ((SJ would 
depend on i.) The two-site Hamiltonian H is stochastic, that is, its non-diagonal 
elements are nonnegative, and the sum of the elements of each of its columns 
vanishes: 

{{s\®{s\)H^Q. (10) 

Here iJ is a 4 x 4 matrix (as the system under consideration has two possible 
states in each site and the interactions are nearest neighbor). The non-diagonal 
elements of H are nonnegative and equal to the interaction rates; that is, the 
element Hp with a ^ (3 \s equal to the rate of change of the state /3 to the 
state a. a and [3, each represent the state of two adjacent sites. For example if 
a = 01 and (3 — 10, then is the rate of particle diffusion to the right. 

The evolution equation of one-point function (n^) (fc-point functions) de- 
pends on both one- and two-point functions (fc — 1-, k- and k + 1-point func- 
tions). Generally this set of evolution equations cannot be solved exactly. If 
one can obtain the state of the system |P) exactly, there is no need to solve this 
set of evolution equations. In [21], a ten-parameter family of reaction-diffusion 
processes was introduced for which the evolution equation of fc-point functions 
contains only k- or less- point functions. We call such systems autonomous. 
The average particle-number in each site has been obtained exactly for these 
models. In [22,23], this has been generalized to multi-species systems and more- 
than- two-site interactions. 

Although generally, one cannot obtain the state of the system |P) exactly, for 
a special choice of initial states and of course with some constraints on reaction 
rates, one may obtain the state of the system |P). There are three families of 



3 



stochastic onc-dimcnsional non-cquihbrium lattice models for which the single- 
shock measures are an invariant subspace of the states of these models [5] . If the 
initial state of these models is a linear superposition of shock measures then one 
can obtain the state of the system |P) exactly. These models are the Asymmetric 
Simple Exclusion Process (ASEP), the Branching- Coalescing Random Walk 
(BCRW) and the Asymmetric Kawasaki-Glauber process (AKGP). 

Let's consider a one-dimensional lattice with L sites. The interaction is 
nearest-neighbor, if the Hamiltonian is of the form 



w = 61 o + 



.i=l 



+ (11) 



where Hi^i^i denotes interaction in the bulk and bi and bL denote the interac- 
tions at the boundaries, 



a —7 J ' \ 5 —(3 

a and 7 (S and f3) are the rates of injection and extraction at the first site (at the 
L'th site). Each site may be occupied or vacant. We represent any configuration 
of the system by the vector \Ea). So the system is spanned by 2^ vectors, \Ea) 
(a = 1, 2, • • • 2^), and any physical state is a linear combination of these vectors 

2^ 2^ 
\^)=Y^^a\Ea), where ^7'a = l. (13) 

a=l a— 1 

Va's are nonnegative real numbers. Va is the probability of finding the system 
in the configuration a. 

It is said that the state of the system is a single-shock at the site k if there 
is a jump in the density at the site k and the state of the system is represented 
by a tensor product of the states at each site as 

where 

It is seen that 

(§|efc) = 1. (16) 

|efc) represents a state for which the occupation probability for the first k sites 
is pi, and the occupation probability for the next L — k sites is p2- The set 
\ek),k = 0, 1, • • • L is not a complete set, but linearly independent. 

There are three families of stochastic one-dimensional non-equilibrium lattice 
models, (ASEP,BCRW,AKGP), for which if the initial state of these models is 
a linear superposition of shock states, at the later times the state of the system 
|P) remains a linear combination of shock state. For these models 

W|efe) = di\ek-i) + d2\ek+i) - {di + d2)ek. (17) 



4 



where di's are some parameters depending on the reaction rates in the bulk, pi 
and p2- So the span of |efc)'s is an invariant subspace of H., the Hamiltonian of 
the above mentioned models. It should be noted that the number of |efe)'s are 
L + I, and any physical state is not necessarily expressible in terms of |e/c)'s. 
For a finite lattice with the injection and extraction at the boundaries besides 
p7|l . there are two other relations. 

n\eo) = Di|ei) -Di|eo) 

n\eL) - D2\eL-i) -D2\eL), (18) 

where D^s are two parameters generally depending on pi and p2, and the reac- 
tion rates. 

Let's assume that the initial state of the system is a linear combination of 
shock state 

L 

|P)(0) = EP'=(0)l^fc)- (19) 

PkS, are not necessarily nonnegative, and so any of them may be greater than 
one. For such an initial state, the system remains in the sub-space spanned by 
shock measures. 

L 

ip)w = Ep'=wiefc)- (20) 

k=0 

Using (|16|l . it is seen that 

L 

fc=0 

But it should be noted that these are not probabilities, pk only expresses the 
contribution of a shock at the site k in the state of the system. Any shock state 
jcfc) can be expanded in terms of \Ea), 

2^ 

|efc) =E^fc-l^->' (22) 

a=l 

where the elements of A are nonnegative. Substituting (|22|l in (|2(J|I and com- 
paring with (|13|) gives 

L 

Va^Y.P'^^k^- (23) 

Here Pa is the probability to find the system in the state \Ea), and so it is a 
nonnegative number. The condition of non-negativeness of probabilities (Vas) 
leads to constraints on pk, see ()23(l . 

The three models are classified as following 

• 1. ASEP- The only non- vanishing rates in the bulk are the rates of dif- 
fusion to the right 0^23 and diffusion to the left W32. For a finite lattice, 
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there may be injection and extraction rates at the boundaries, a, and 7 
are the injection and extraction rates at the left boundary, and 6, and (3 
the injection and extraction rates at the right boundary. In this case the 
densities can take any value between and 1 (pi ^ 0, 1 p2 7^ 0, 1). rfi, 
and d2 are 

dl = [W23 - W32) 

P2 - Pi 

, P2(l - P2) I . 

di = (W23-W32). (24) 

It should be noted that the densities pi, and p2 are also related through 

P2(l - pi) _ W23 



Pl(l-P2) W32 

So 



(25) 



d\ = —W23, d2 = —C032- (26) 

P2 Pi 

The rates of injection and extraction at the boundaries are also related to 
the densities pi, and p2- 

pi(l - pi)(a;23 - W32) = a{l - pi) - 7pi, 
-P2(l - P2)(W23 - W32) = ^(1 - P2) - (3p2- (27) 

The parameters Di, and D2, are obtained to be 

-Dl = 02 H = ^2 H W23, 

Pi - P2 Pi 

D2 = dl- ^^^-P^)-^P^ =dl + ^-.32. (28) 
Pi - P2 P2 

2. BCRW - The non- vanishing rates arc coalescence (CJ34, and UJ24), Branch- 
ing (a;42, and (JJ43) and diffusion to the left and right (a;32, and a;23)- The 
injection rate at the right boundary 6 should be zero. The density pi can 
take any value between and 1, but p2 should be zero. These parameters 
are related through 

W23 1 - Pi 



W43 Pi 

W23 _ '<^24 + W34 

W43 a;42 + a;43 



(29) 



The rates of injection and extraction at the boundaries are also related to 
the densities pi, and the reaction rates, 

P1W23 - Pi(l - Pi)w32 - P?W34 = a(l - pi) - 7pi. (30) 
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The parameters di, and d2 are 




(31) 



^1 = — , 

Pi 

D2 = dl - d2{l - pi) + 0. 



(32) 



• 3. AKGP - The non- vanishing rates are Death (0^12, and cui^) and Branch- 
ing to the left and right {0142, and 0)43), and also diffusion to the left, pi 
should be equal to one, and p2 should be zero. The extraction rate at the 
left boundary 7, and the injection rate at the right boundary S should be 
zero. The hoping parameters are rfi = wis, d2 = ^43, and finally Di = a, 
and D2 = /?. 

Interchanging pi and p2, is nothing but exchanging left and right. Changing pi 
to 1 — Pi is particle hole exchange. 

3 single-shock 

3.1 Shock on an infinite lattice 

Here, we want to consider the evolution of shock measures on an infinite lattice. 

Here je;) stands for a state with a shock at the site I. If the initial state is 
a linear combination of shock measures, at later times the state of the system 
should be also expressible in terms of shock measures. 



Calculating p^, one can obtain any correlation functions of number operators. 




(33) 



k=—oo 



00 




(34) 



where 




(35) 



So 



{ni) = P2 + [pi 



P2)Bi{t) 



(36) 



where 



00 




(37) 



k=i 
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All other correlation functions of number operators can be obtained in terms of 



{rii) = P2 + [pi - P2)Bi{t) 
(mnj) = pI + {pi - p2){p2Bi(t) + piBj{t)} 
{niUjUk) = pI + {pi - p2){plBi(t) + p2PiBj{t) + p\Bk{t)) 



where i < j < k < ■ ■ ■ . 

The evolution equation for the system is 



(38) 



n\F) = ^|P). (39) 

Knowing the action of Ti. on |efc), one can obtain the evolution equation for 
Pk{t)- 

Pk = dipk+i + d2Pk-i - (di + d2)pk- (40) 

Here we have used the linear independence of |efc)'s. Let's define the generating 
function 

oo 

G(z,i) := ^'W^'- (41) 

/ — — oo 

Then the evolution equation for G{z, t) is 

G^[diz-^ + d2Z-{di + d2)]G, (42) 
the solution for which is 

G(z, t) = e^d^z-^+d^--(di+d2)]tQ^^^ 0). (43) 

G(z, 0) can be determined using (|41|l and contributions of the shock measures 
in the initial state. The coefficients for the Laurent expansion of the generating 
function are pk{t)s; 

°° \ (fc-™)/2 

p,(f) = e-('^^+'^^)* ^ if] Ik-m{2Vd^2t)pm{Q). (44) 

m— — oo ^ ^ ' 

The above result is first obtained in [4]. At large times, 

»<"nS Ti ■ <^^' 

It is seen from the above equation that if ^2 < di , the contribution of the shocks 
at the rightmost sites tend rapidly to their final value, and obviously for d2 > di 
the contribution of the shocks at the leftmost sites arrive earlier to their final 
values. This expression seems to be unbounded for k — > ±00. For any fixed t, 
this is true. However, it simply means that in order that this term represents 
the leading term for some fc, t must be greater than some T, which does depend 
on k. 
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3.2 Shocks on a lattice with the boundary 

Now let's first consider a one-dimensional lattice with L sites. There are in- 
jection and extraction at the boundaries. Because of the boundary terms p7|l 
changes to 

Ti\ek) = di|efe_i) -I- d2|efc+i> - (di -f d2)|efc) 
H|eo) = A|ei)-A|eo) 

Wl) = D2\eL-i) - D2\eL) (46) 

The evolution equation for p^s in the bulk are the same as that of infinite lattice, 
but here one should take care of boundary terms. 

Pk = dipk+i + d2Pk-i - {di + d2)pk, k ^0,l,L - l,L 

PL-l = D2PL + d2PL-2 ~ (dl + d2)pL-l, 

Pi = dip2 + DiPq - {di + d2)pi, 

PL = d2PL-l- D2PL, 

Po = dipi - DiPq. (47) 

Let's first consider the stationary case. These set of equations can be solved 
easily, 

Pk-i^r\^y-H^r-"Pi (48) 

di D2 Di 

For any finite L, all the p^s can be obtained in term of pi, and pi is also obtained 
using the normalization condition (|21|l . In the thermodynamic limit (L — )■ 00), 
it is seen that for c?2 > ^i, finiteness of pk for large k, leads to vanishingly 
small pi . But as di exceeds ^2 , Pk for large k becomes vanishingly small in the 
thermodynamic limit. This is the static phase transition previously mentioned. 
The static phase transition is controlled by the reaction rates in the bulk and 
is independent of the rates at the boundaries. It is a discontinuous change of 
the behavior of the derivative of the stationary value of pk at the end points, 
with respect to the reaction rates. This phase transition is a well known first 
order phase transition for the familiar ASEP, and also its extensions [20]. The 
BCRW, and AKGP models shows the same phase transitions. In the stationary 
state, the system can be found in the low-density or in the high-density phase, 
or in the other words the contribution of leftmost or rightmost shocks may be 
negligible. 

Now, return to the dynamics of the system. To proceed let's first introduce 
a change of variable 

qk ■■=Pk, k ^0,L 
qo := {Di/d2)po, 

qL ■■= {D2/di)pL. (49) 
Using this change of variable H47|) recasts to 

qk = diqk+i + d2qk-i - [di + d2)gfe, fc 0, L 
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<1L = 

90 = 



diDi 



■QL-i - D2qL, 
qi - Diqo. 



(50) 



To find the relaxation of the system towards its stationary state, one should find 
the greatest nonzero eigenvalue of the operator h, defined through =■ h^^^qi- 
The eigenvalues and eigenvectors of h have been denoted by E and C e , respec- 
tively. Expanding the vector q in terms of C^'s, and regarding the completeness 
and linear independency of C^'s, one arrives at 



E Ck = diCfe+i + daCfe-i - (di + d2)Ck 

d2D2 

E Cl = — ; — Cl-1 — D2CL, 
E Co — — ^Ci — DiCq. 



(51) 



The solution to the above equations is 
where Zi's satisfy 



diZi H . 

z,. 



E^-{di+d2 

Then Z1Z2 = d2/di. The second and third equations of (|51|l take the form 



(52) 
(53) 



{E + D2){az^ + bz^) 



d-,D 



di 



hZ' 



L-l\ 



0, 



{E + D^){a + 5) - ^^{azi + 6Z2) = 0. 

02 



(54) 



or 



iE + D2)zi 



d2D2 



di 



{E+D2)Z2- 



d2D2 



di 



L-l 



6 = 0, 



diDi 
E + Di- -\^zi 



diDi 
E + Di- ^r^Z2 



6 = 0. 



(55) 



To have nonzero solutions for g^'s, these equations should have nontrivial solu- 
tions for a, and 6, which means that the determinant of the coefficients should 
be zero, 



E + D.-'^-l^z, 



E + D,-^Z2 



{E + D2)Z2 " 

{E + D2)zi 



do Do 



di 

d2D2 



L-l 



di 



(56) 
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Performing the change of variable Zi =: \fd2fd\ Xi, leads to 



E = -{di + da) + Vd^{x^ + 
and X1X2 = 1- The equation H56|l changes to 



E + Di - Dixiy%jd2 {E + D2) - D2Xiy^d2/di 



-2L+2 



{E + Di)xi~ Diy/di/d2 iE + D2)xi~D2y/d2/di 



(57) 



(58) 



Equation (|58|l can be written as a polynomial equation of order 2L, so it has 2L 
solutions. Two obvious solutions of the are x = ±1. But, these generally do 
not correspond to eigenvalues and eigenvectors. In fact for these solutions, xi 
and X2 are the same, so that should be modified to Ck = {a + bk){—l)^ and 
it is not difficult to see that these do not fulfill the boundary conditions unless 
a — b — 0. It will be shown that two cases may occur, either both solutions are 
phases then |a:i| = \x2\ = 1, or both of them are real . Except for the solutions 
±1, one of the real solutions is greater than one, which we take the be Xi. For 
the phase solution, x = exp(i0) 



E = -{di + d2) + 2y^did2 cos 6*. 



(59) 



Among these set of eigenvalues the biggest one is £' = —{^/di — v^)^- So if 
there is no other solution except for the phase solutions, the relaxation time is 



1 



To 



(VHT 



(60) 



Now, let's search for the real solutions, if those exist. Of course, by real solutions 
it is meant real solutions besides the trivial solutions ±1. If a; is a solution to 
H58(l . x~^ is another solution to it. So it is sufficient to seek the solutions with 
I a; I > 1. In the thermodynamic limit, and for |a:;| > 1, H58|) is simplified to 



{E + D2)x-D2y/d2/di {E + Di)x - Diy/d^ 



(61) 



If such a solution exists, then 



-{di 



^/d^{x + x-^). (62) 

for any positive x, E obtained from H62|l is greater than the eigenvalue obtained 
from H59|l . and so the system relaxes to its stationary state slower. If all of the 
solutions for (|58|) are phases we call the system is in fast phase, and if there 
is real solution, the system is in the slow phase. Depending on parameters, 
there may be more than one real solutions, and the system may be in the slow 
or slower phase. The solutions of first bracket in l|?>T|l is cc = y/d2/di, and 
(rfi — D2) / y/dida- The first solution gives E = Q, which is related to stationary 
state. We obtained these solutions assuming > 1. So, x = {di — D2)l \Jd\d2 
is a solution provided that 

dx 



D2 



Vdid^ 



> 1 



(63) 
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Similarly there may be a solution for the second bracket in Utilf) . x = {d2 — 
Di)/y/did2, provided that 

B:='-^>1 (64) 

The above results are summarized in the figure 1. Depending on the magnitudes 
of A, B, there are five regions. If both A, B are less than one, then the relaxation 
time is given by (|60() (region /). Changing continuously the reaction rates at the 
bulk or at the boundaries or the densities, the relaxation time may changes 
discontinuously. This is the dynamical phase transition. If one of A^ and B 
become greater than one, then the relaxation time is given by T4 or re 

di - D2 

D2{di-D2-d2) 

^« ^= D,id2~D,-d,y 

If ^ < 1, and B > 1 the system is in the phase //. If ^ > 1, and B < 1 the 
system is in the phase ///. And finally if it is possible that both A, and B 
become greater than one, the relaxation time will be given by max(T4, re). In 
the phase IV, the relaxation time is re, and in the phase V, the relaxation time 

is TA- 

Let's study each of the three models ASEP, BCRW, and AKGP separately. 

• 1. ASEP - To study the phase structure of this model, we should obtain 
A, and B. 

Let's assume ^32 < ^23- Then LU32 < ^/oJ23i^, and obviously 

S 



W32 < ^/ 0^231^32 ■ 

Now using the second equation of H28(l . one arrives at 



dl - D2 < \J^2Z^Z2 Vdid2, ^ < 1. 

Now, let's assume a;23 < a;32. Using (|27|l . and the fact that (3 is the 
extraction rate, and should be positive, one gets 

S 

^32 < ^23- 

P2 

But W23 < ^^^23^^32, so 

di - D2 < ^^23^32 = \/dld2, =^ ^ < 1. 

With the similar reasoning, one gets B < 1. Therefore the only region in 
the space of parameters A, and B available for the ASEP is the region I. 
So, there is no dynamical phase transition for ASEP, and the relaxation 
time is given by 

- . (66) 



(Pl V^23 - P2y/i032y 
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• 2. BCRW- Let's first define 



n := UJ32 + (67) 
i - pi 

Using H3U|I . and tlie fact tfiat 7 is tlie extraction rate, and sliould be 
positive, one gets 

^iinZll > _ (1 _ p^)n. (68) 
Pi 

Ttiree cases may occur. 



— UJ23 < fi, and a(l — pi)/pi < W23 ~ (1 " Pi)aA^23^^- Tlren ^ < 1, and 
B > 1. Tlierefore, tlie system is in tlie phase //, and the relaxation 
time is given by rg. 

— W23 > n, and P < a;23 - (1 - Pi)^- Then ^ > 1, and B < 1. 
Therefore, the system is in the phase ///, and the relaxation time is 
given by r^. 

— Otherwise the system is in the phase /, and the relaxation time is 
given by tq, (jHOJl. 

The regions with ^ > 1, and S > 1 is not available for the BCRW. 
In summary, for the model BCRW, the parameter space is restricted 
to the regions /, //, and ///. So, there are three distinct phases 
available for the system, and this model may experience dynamical 
phase transitions. 

• 3. AKGP- For this model we have 




(69) 

W13W43 

If Wi3 > W43, Then B < 1. Depending on the rates UJ13 and W43, and /3, the 
system may be in the phases /, or ///. And if CJ13 < ^43, Then ^ < 1. 
Depending on the rates ^13 and 0-143, and a, the system may be in the 
phases /, or //. For the model AKGP, the parameter space is restricted 
to the regions /, //, and ///. So, there are three distinct phases available 
for the system, and the model AKGP may experience dynamical phase 
transitions. 



Double shock 

define the state of a double shock on a one-dimensional lattice with L sites 
|e„,fe) m + k<L (70) 
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where 

Pi J \ P2 J \ P3 



(71) 



|em,fc) represents a state for which the occupation probabihty for the first m sites 
is pi, the occupation probabihty for the next k sites is p2, and the occupation 
probabihty for remaining sites is ps. It should be noted that this state represents 
a double shock, one shock at the site m, and the other one at the site m + k. We 
call k the width of double-shock. We also assume that the three densities pi, 
P2, and p3 are different. The evolution of the shocks in ASEP starting from a 
measure with several shocks and extra particles at the shock positions is studied 
in [4]. In that model, the measure is 

\em.k) = u'^"" ^ Z ^ v'^'' (E) Z (E) w'^^^-''-"'\ m + k<L, (72) 

where 

It is shown that with a special choice of densities, time evolution equation of 
this model is similar to that of some random walkers on a lattice. 

Although the states \em,k), {m, = 0, 1, • • • ) do not construct a complete 
basis, they are linearly independent provided that k the width of double-shock 
can be zero only when m = or m = L. This means that there should be 
at least one site with the occupation probability p2 between the sites with the 
occupation probability pi, and the sites with the occupation probability ps. For 
the single-shock linear independency of \ek)s are obvious, but here m, v, and 
w are not linearly independent. So it is not obvious that |em,fc)'s are linearly 
independent. To prove |em,fc)'s are linearly independent, one must show that if 

L L-k 

J2 = 0, (74) 

fc=0 m=0 

then all a'^'^'s must be zero. Let's define 

u := [pi pi~l) , V := [p2 P2-1) , w := (ps ps - l) ■ (75) 
Multiplying equation (|74|l from left hand side by js)**^, one arrives at 

EE«'"' = 0- (76) 

fe=0 m=0 

Prime on the first summation denotes that fc = 0, only when m = 0, or m = 
L. Now, multiplying equation 174|l from left hand side by (sj*^^^^ E) u, and 
w (E) {sl^'"^^ one arrives at 

L L—k 
k=0 m=0 
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L L-k 

J2 J2 - = ^- c^^) 

fe=Om=0 

(O together with lO gives = a"" = 0. Now, define Z^-" through 

/™ := w®'' ® s®" ® (i-'--") ^ (78) 

where r, and n are non-negative numbers and n + r < L. Multiplying equation 
ifT^ from left hand side by f^'^ and using f'^'^\em,k) = ^m^l (fo'^ ^ 7^ 0), one 
arrives at a*"^ = 0. Similarly, it can be shown that all the coefficients a™'' in 
ifTHl should be zero. So, the set |em_fc)'s are linearly independent. 

It was previously mentioned that there are three families of stochastic one- 
dimensional non-equilibrium lattice models, (ASEP,BCRW,AKGP), for which 
if the initial state is a linear superposition of single-shock measures, at later 
times the state of the system |P) remains a linear combination of shock mea- 
sures. Among these, ASEP is the only model for which double-shocks forms an 
invariant subspace, which means that if the initial state is a linear superposition 
of the double-shock measures, at later times the state of the system |P) remains 
a linear combination of double-shocks. It can be shown that for the case of 
ASEP and on a one-dimensional lattice with infinite sites, 

Ti-lem^k) = rfl|em-l,fc+l> + d2\e,n+l,k-l) + d3\e,n,k-l) + rf4|em,fc+l) 

-idi + d2 + d3 + di)\e^,k), k>2 

Ti-lem,!) = di\em.-i,2) + d4\em,2) ~ [di + d4:)\em,i) , (79) 
where 



and 



^23 _ P2(l - Pi) _ P3(l - P2) 
W32 Pl(l-/32) P2(l-P3)' 



(80) 



di 


= (W23 - 


- t^32) 


d2 


= (W23 - 


~ t^32) 


ds 


= (W23 - 


- W32) 


C?4 


= (W23 - 


~ ^^32) 



P2 - 


Pi 


P2{1- 


- P2) 


P2 - 


Pi 


P2(l- 


~ P2) 


P3 - 


P2 


Psil' 


- Pi) 


P3 - 


P2 



(81) 

If the initial state is a linear combination of double-shocks, then 

00 00 

l^)W= I^P^feWlemfe), (82) 

jri=— 00 k=l 

where pm,fc is the contribution of the double-shock, rak in the state of the system. 
Using . (182(1 . and also the linear independency of |em,fc)'s, one can obtain the 
evolution equation for pmk's. 

Prn.k — diPrn+l,k-l + '^2Pm-l,fc+l + C^3Pm,fc+l + C?4Pm,fe-l 
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-(di +d2+d3+ d4)p,n,k, 



k>2 



Pra,l = d2Pm-l,2 + " (^1 + dAjPm,! (83) 

Let's define 

oo 

qk-= ^ Pmk- (84) 

m— — oo 

Then qk is the contribution of all double-shocks with the width, the distance 
between two shocks, k in the state of the system. Conservation of the probability, 
leads to 

oo 

(85) 

k=l 

The evolution equations for q^'s are 

qk = DiQk-i + D2qk+i ~ {Di + £'2)%, k>2 

qi = D2q2-Diqi, (86) 

where 

Di-~di+di, D2:=d2 + d3. (87) 
Di, and D2 can be written in terms of pi and p2, and the diffusion rates. 

Di = (t^32-W23){ 1 } 

Pi ~ P2 P2- P3 

D2 = (-3.--23){^^^i^ + ^^^i^} (88) 
Pi ~ P2 P2- P3 

Subtracting these 

D - D ^ P2(l-P2)(l + f)(l-f-)^ 

' ' C^23[(1-P2)+FP2][^^(1-P2)+P2] ^ ' 

where F :— 1^32/^^23, it is seen that Di < D2- One can easily obtain the steady 
state solution for g^'s. 

,1 = (90) 

To obtain the second equation of (|90|l . we have used H85() . Using the fact that 
Di < D2, it is seen that (7^. goes to zero for large fc's. This means that at the 
stationary state the contributions of the double-shocks with larger width are less 
and the main contribution comes from thin double-shocks. This is reminiscent 
of double-shocks in Burgers equation. In fact, in the hydrodynamic limit double 
shock is not stable and converges to single shock. Here, we show that in the 
microscopic level, and in the thermodynamic limit, the stationary value of the 
contribution of double shocks with larger width become vanishingly small. 
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Figure 1- Exact phase diagram of the ASEP, BCRW, and AKGP, in the A, 
B plane. For the ASEP, the parameter space is restricted to the region /. For 
the models BCRW, and AKGP, the parameter space is restricted to the regions 
J, II, and ///. 

Figure 1- Exact phase diagram of the ASEP, BCRW, and AKCP, in the A, 
B plane. For the ASEP, the parameter space is restricted to the region I. For 
the models BCRW, and AKGP, the parameter space is restricted to the regions 
J, //, and ///. 
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